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ABSTRACT 

We present results from an intensive VERITAS monitoring campaign of the 
high-frequency peaked BL Lac object lES 1218+304 in 2008/2009. Although 
lES 1218+304 was detected previously by MAGIC and VERITAS at a persistent 
level of ~6% of the Crab Nebula flux, the new VERITAS data reveal a prominent 
flare reaching ~20% of the Crab. While very high energy (VHE) flares are quite 
common in many nearby blazars, the case of lES 1218+304 (redshift z — 0.182) is 
particularly interesting since it belongs to a group of blazars that exhibit unusu- 
ally hard VHE spectra considering their redshifts. When correcting the measured 
spectra for absorption by the extragalactic background light, lES 1218+304 and 
a number of other blazars are found to have differential photon indices F < 1.5. 
The difficulty in modeling these hard spectral energy distributions in blazar jets 
has led to a range of theoretical 7-ray emission scenarios, one of which is strongly 
constrained by these new VERITAS observations. We consider the imphcations 
of the observed light curve of lES 1218+304, which shows day-scale ffux varia- 
tions, for shock acceleration scenarios in relativistic jets, and in particular for the 
viability of kiloparsec-scale jet emission scenarios. 

Suhject headings: BL Lacertae objects: individual (lES 1218+304 = VER J1221+301) 
— galaxies: active — gamma rays: galaxies 
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INTRODUCTION 



The very high energy (VHE) 7-ray sky has seen a dramatic increase in the number 
of sources, with 26 detected active galactic nuclejl] (AGNs) over redshifts of 0.0018 
-0.536. All, except two, are blazars with their jets closely aligned to our line of sight. 
Blazars as a class are co mprised of BL Lacertae objects, optically violent variables and 



other radio-loud quasars (lUrry fc Padovani 



19951 ). Non-thermal emission dominates the 



broad-band continuum spectra of blazars. 

The characteristic double-peaked spectral energy distribution of blazars is usually 
attributed to synchrotron radiation in the radio to X-ray wavebands, and inverse- 



Compton up-s cattering of so 



TeV energies (IMaraschi et al. 



t target photons by relativistic electrons at GeV and 



1992 



Bloom fc Marscher 



19961 ). Alternative emission 
scenarios involve hadronic in teractions producing neutral pions which decay into p hotons 



( iMannheim fc Bierman 



I992I ) and proton synchrotron radiation (lAharonianl |2000| ) . The 
latter models provide a link between AGNs and possible sources of ultra-high-energy 
cosmic rays. A better understanding of the underlying acceleration and emission processes 
requires a larger sample of sources, which is now being provided by the Fermi Gamma-ray 
Space Telescope (Fermi) at 0.1 GeV to hundreds of GeV and by the imaging atmospheric 
Cherenkov telescopes (lACT), namely the MAGIC, H.E.S.S. and VERITAS Observatories, 
that provide coverage primarily in the 100 GeV — 10 TeV regime. 

To date, blazars detected at VHE energies are predominantly high-frequency peaked 



BL Lacertae objects (HBLs). Nevertheless, a few objects with 



low-frequency-peaked BL Lacertae (LBLs; (lAlbert et al. 



2007 



ower peak energies, 



Teshima et al. 



intermediate-frequency-peaked BL Lacertae objects (IBLs) (see 



Acciari et al. 



namely 



2008V) and 



2008a 



^For updates see |http: //tevcat. uchicago.edu 



- 6 - 



2009bl : lOng et al.l |2009| ) were detected recently by lACTs at TeV energies. The 7-ray 
luminosity of LBLs and IBLs generally peaks at sub-GeV to tens of GeV in energy, and 
they are therefore easily detectable by Fermi while their detectability in the TeV regime by 
lACTs is more difficult due to their soft VHE spectra. 

Among the HBLs, which are typically the domain of the TeV telescopes, a number 
of objects exhibit unusually hard intrinsic power law energy spectra (dN/dE ~ E"'^') 
after correcting for the 77 e"*" e~ absorption by the cosmological diffuse EEL 
radiation field. While the measured spectral indices, Fm, of these blazars (lES 1101-232, 
lES 0347-121, lES 0229+200, lES 1218+304, RGB J0710+591) range from 2.5 to 3.1 



see 



Aharonian et al. 



2006ai l2007al Jb 



Albert et al. 



2006 



Acciari et al. 



2009a 



Perkins et al. 



20091 ). the absorption-corrected spectral in dices suggest very hard intrinsic source spectra 



in the VHE regime with T; < 1.28 ± .28 ( 



limits to the EEL from galaxy counts f lLevenson fc Wright 



Krennrich et a. 



20081 ). 



20081) . based on recent lower 



Similar hard energy spectra are found by Fermi, where, for example, the measured 
spectral indices for RGB J0710+591 and lES 1440+122 are found to be T = 1.21 ± 0.25 



2009ah . While the 



and 1.18 ± 0.27 over energies from 0.2 GeV to tens of GeV (lAbdo et al.l 
findings of very hard VHE spectra are based on EBL constraints from galaxy counts, the 
Fermi spectra directly resemble the intrinsic source spectra since absorption by the EBL is 
minimal for the covered energy range. The photon spectrum of lES 1218+304 measured 
by Fermi yields a spectral index of F = 1.63 ± 0.12 betwee n 0.2 GeV and ~3 00 GeV, but 



absorption may already play a role at the high-energy end (lAbdo et al. 



2009ah 



Malkov fc Drury 



2OOII ) generally 



Diffusive shock-acceleration theory (for a review see 
yields a limit to the spectral index of GeV - TeV photon spectra resulting from inverse 
Compton scattering of Fi > 1.5. Only recently, numerical studies by Stecker et al. (2007) 
indicate that sufficiently hard electron spectra could be generated by diffusive shock 
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accelerati on in relativisti c shoc ks, for the production of photon spectra 1.0 < Fi < 1.5 
However 



Bottcher et al. 



(120081 ) suggest that even for a hard-spectrum electron population, 
in the framework of a synchrotron self-Compton (SSC) scenario, the resulting GeV-TeV 
7-ray spectra should experience substantial softening from Klein-Nishina effects making 
7-ray spectra with Fi < 1.5 less likely. 

Other approaches to explain the hard 7-ray spectra are offered by ad hoc assumptions 
about the electron distribution, an additional absorption component in the source, or a 
postul ate of new physics describ ing the propagation of 7-rays on extragalactic distance 



scales. 



Katarzyhski et al 



(I2OO6I ) invoke a high low-energy cutoff in the electron distribution 



that could give t he ap pearance of a hard 7-ray spectrum for a given energy regime. 



Aharonian et al. 



( I2OO8I ) show that 77 absorption in the source due to a narrow-band 



emission component f rom t he AGN could lead to unusually hard VHE spectra. Finally, 



Sanchez-Conde et al. 



(120091 ) suggested an axion like particle (ALP) that would distort the 
7-ray spectrum through ALP/photon mixing on cosmological distances in the presence of 
intergalactic magnetic fields. 

A more easily te stable model that involves known physics was recently proposed by 



Bottcher et al. 



(120081 ). The authors attribute the hard photon spectra to inverse-Compton 
up-scattering of ambient photons from the cosmic microwave background (CMB), occurring 
in a kiloparsec-scale jet. In this case, a substantial fraction of the jet power is transported by 
hadrons to the outer regions of the jet, where it is dissipated into ultra-relativistic electrons. 
Inverse-Compton scattering in the Klein-Nishina regime is no longer the limitation at the 
highest energies since low-energy target photons from the CMB are abundant. The low 
magnetic field in the large-scale jet also avoids the overproduction of synchrotron radiation 
and allows Compton emission to dominate. Furthermore, the synchrotron emission in 
the radio to X-ray wavebands is assumed to originate on sub-parsec scales. A similar 
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mechanism was previously suggested to explain the hard X -ray spectra observed from the 
kiloparsec-scale jets of radio quasars ( iTavecchio et al.l 120071 ). 



To date, all of the blazars exhibiting very hard spectra appear to emit at a baseline 
level that is consistent with this picture. However, it is also possible that the sensitivity 
of the current lACTs is the limiting factor in detecting day scale variations and that the 
underlying emission may contain flares. The question as to whether or not the 7-ray 
emission from lES 1218+304 consists of flares or corresponds to a constant baseline emission 
level was one of the main motivations for VERITAS to monitor this object over a period of 
5 months in the 2008/2009 season. The other motivation for deep exposures of relatively 
large redshift {z = 0.1 — 0.2) blazars at TeV energies is to provide better constraints on the 
EBL spectrum. The observed VHE photon spectrum can be used to constrain the EBL in 
the near to mid-IR. This is particularly promising for hard-spectrum blazars that provide 
the best sensitivity to any possible absorption feature. 



2. OBSERVATIONS AND DATA ANALYSIS 

VERITAS is an array of four 12 m lACTs located at the base camp of the F.L. 
Whipple Observatory in southern Arizona (31?68N, 110?95W, 1.3 km above sea level). Each 
VERITAS focal plane instrument has a 3?5 field of view and consists of 499 photomultiplier 
tubes. The stereoscopic system has been fully operational since mid-2007 and is capable 
of detecting sources with a fiux of 1% of the Crab Nebula in < 50 hr of observation. In 
Summer 2009, Telescope 1 was relocated and the improved array configuration allows us to 
detect a point source at 1% of the Crab Nebula fiux under 30 hr of observations. The large 
effective area (~ 10^ m^) of the array enables VERITAS to be sensitive over a wide range 
of energies (from 100 GeV to 30 TeV) with an ene rgy resolution of 15 % — 20% above 300 



GeV. For further details about VERITAS, see e.g.. 



Holder et al 



fl2006h . 
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lES 1218+304 was observed by VERITAS from 2008 December to 2009 May. All 
data were taken in wobble mode where the source is positioned 0?5 from the center of the 
camera in order to allow simultaneous background estimation. A quality selection was 
applied to the data to reject observations affected by poor weather conditions and other 
anomalies, resulting in a total live time of 27.2 hr for analysis. The mean zenith angle of 
the observations was 12°. 

The data are analyzed using standard VERITAS calibration and analysis tools. 
Initially, the observed data set is calibrated by flat-fielding the camera gains with nightly 
laser runs. Following cali bration and c leaning, the shower images are parameterized using a 



second-moment analysis (IHillas 



19851 ). Images from individual telescopes are combined to 
reconstruct both arrival direction and impact position for each event. The 7-ray-like events 
are separated from the cosmic ray background events using cuts on ihe mean-r educed- scaled 



Aharonian et al 



2006b for a full 



length and mean-reduced- scaled width parameters (see 
description of the parameters). The cuts for 7-ray/hadron separation are optimized a priori 
for a source with a 10% Crab-like fiux. A circular region of radius 9 is applied around 
the target to estimate the ON-source eve nts. The backgrou nd or OFF-source counts are 



estimated with off-source reflected regions (IBerge et al. 



20071 ) . The sig nifican ce of any excess 



(ON-OFF) events is calculated following Equation (17) in Li & Ma ( 1l983l ). The angular 
distribution of excess events from the direction of lES 1218+304 (Figured]) is comparable 
with the expectation from the Monte Carlo (MC) simulation of a point source. 

MC simulations are used to generate a lookup table for estimating the energies of the 
primary 7-rays as a function of the image size and impact parameters. For spectral analysis, 
a looser cut {spectrum cuts) on the arrival direction of 7-rays (a region with 9"^ = 0?03) 
is applied along with three off-source reflected-regions. This has the benefit of retaining 
a higher number of 7-ray candidates in addition to reducing the systematic uncertainties 
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associated with the MC detector model. Both the standard cuts (6*^ = 0?0225) and the 
spectrum cuts yield compatible results. Consistent results are attained with independent 
analysis packages. 



3. RESULTS 

Table [T] summarizes the results of the VERITAS observations of lES 1218+304. For the 
spectral analysis, we report an excess of 1155 events with a statistical significance of 21.8 
standard deviations, a, from the direction of lES 1218+304 during the 2008-2009 campaign 
(2808 signal events, 4959 background events with a normalization of 0.33). Figure [2] shows 
the corresponding time-averaged differential energy spectrum. The spectrum extends from 
200 GeV to 1.8 TeV and is well described (x^/dof = 8.2/7) by a power law. 



— = (11.5 ± 0.7) X 10-12 ( ^ )-3.o7±o.o9.... cm-^s-^TeV-i (1) 
dE ^ ' ^500GeV^ ^ ' 

The errors bars are statistical only. 

Figure E] shows the light curve of lES 1218+304 for the 2008-2009 season (binned 
by night). The average integral photon flux above 200 GeV is $(E > 200 GeV) = 
(18.4 ± 0.9 Stat) X 10^12 cm^^ s^^ corresponding to ~7% of the Crab Nebula flux above the 
same threshold. The light curve from lES 1218+304 shows strong nightly flux variations 
between 2009 January 25 and 2009 February 5. The highest 7-ray flux is observed on MJD 
54861 (2009 January 30) and is ~20% of the Crab Nebula flux. The inset within Figure |3] 
shows the light curve above 200 GeV for nights with the highest flux. A fit of a constant 
to the emission during the entire season yields a x^/dof of 106/21 (Figure [31 dashed dotted 
line), corresponding to a very small probability (~ 10^^^). The data taken during the 
active period in January/February of 2009 indicate variability on a time-scale of days. 
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In order to investigate possible changes to the spectral shape during the flaring activity, 
the total data set is divided into high-state and low-state samples. The two nights with 
the highest flux during flaring activity (January 30 and 31) are considered the high state, 
while the remaining nights are considered the low state. Both the high-state and low-state 
spectra are well described by a power law, and Table |2] summarizes the spectral results. 
We do not see any evidence for a change in the spectral index during the increase in the 
absolute flux level from lES 1218+304. 



DISCUSSION AND CONCLUSIONS 



A prominent feature in the emission properties of blazars is the observed variability on 
time scales ranging from minutes to years. In 2006, Suzaku detect ed a prominent X -ray flare 



from lES 1218+304 revealing a hard-lag variation in the X-rays (iSato et al. 



lES 1 218+304 was flagged as variable in the Fermi bright so urce catalog 



2008). While 



(lAbdo et al 



2009aj). The 



2009bl ). no variability was detected in the Fermi light curve (lAbdo et al. 
increased flux measured by VERITAS is the flrst signiflcant observation of flaring activity in 
the VHE emission from lES 1218+304. Among the other hard-spectrum VHE blazars (lES 
1101-232, lES 0347-121, lES 0229+200, 1ES1218+304, RGB J0710+591) measured to 
date, lES 1218+304 is the flrst clear example of signiflcant VHE variation from a baseline 



flux. Consequently, the strong flare will challenge existing spectral mode 
that assume a constant VHE emission from the source (see, for example. 



s of lES 1218+304 



Riiger et al. 



20091 ). 



It is generally assumed that the VHE emission from blazars is produced in the 
relativistic jet where the highly collimated and Doppler-boosted emission is responsible for 
the high luminosities and short variability timescales. The size of the emission region, R 
is constrained by the causality argument to R < ctyar^/i^ + z) where 6 is the relativistic 
Doppler factor and t^ar is the observed 7-ray variability timescale. An estimate of the 
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flux-doubling time {t^ar ^ 1 day, see Figure [3]) limits the size of the emission region to 
R6-^ < 2.1 9 X 10^5 m^ 0.71 x lO^^ pc. The Doppler factor {6 ~ 20) typically derived 



for blazars ( iMarscher 



20061 ) implies that R < 0.01 pc. 



While the hard emission spectrum during the quiet state is consistent with the CMB- 
inverse-Compton interpretation, the requirement for a large emission region (~pc scale) 
predicted by the model is excluded by the R < 0.01 pc constraint imposed by the variability. 
Therefore, the flaring behavior from lES 1218+304 implies that the CMB-inverse-Compton 
model of emission in the kiloparsec-scale extended jet is unlikely to be the sole explanation 
for the extreme hardness in the intrinsic spectrum of lES 1218+304. 

Although our observations of lES 1218+304 are consistent with a baseline level of 
VHE emission that could have its origin in an extended emission region, the variability 
described here indicates that comparable flux with a comparably hard spectrum is emitted 
from the compact regions of this source. Further observations of lES 1218+304 are required 
to determine if the hard-spectrum emission can be explained as originating from the 
kiloparsec extended jet via the CMB-inverse-Compton model. Consequently, long-term 
VHE monitoring with ground-based 7-ray telescopes may be crucial to unraveling the 
emission mechanism in distant the hard-spectrum VHE blazars. 
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Fig. 1. — Distribution of 6*^ for excess events ON-OFF from the observation of lES 1218+304. 

The dashed-dotted hne shows the boundary of the region for the spectrum cuts (9^ = 0.03°). 
The sohd curve indicates the expected 9"^ distribution from a point source and provides a 
good fit to the data. 
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Fig. 2. — Observed differential energy spectrum of lES 1218+304 from tlie 2008-2009 data 
set. Tlie solid line shows the power law fit to the data (filled circles) with a spectral index 
value of —3.07 ± 0.09 Stat- This agrees, within statistical uncertainties , with the previous 



VERITAS measurement of —3.08 ± 0.34 stat (dashed line) ( lAcciari et al. 



2009ah . 
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Fig. 3. — Night-by-night VHE hght curve for lES 1218+304 is shown as measured from 
2008 December 29 (MJD 54829) to 2009 April 23 (MJD 54944). The open circles represent 
the integral flux above 200 GeV, 0(E > 200 GeV) from lES 1218+304 assuming a spectral 
shape, dN/dE oc with F = 3.07. The dashed-dotted line corresponds to the average 
integral flux of $(E > 200 GeV) = 18.4 ± 0.9stat x 10"^^ ^^^-2 r^^^ j^g^^ ^^^^^ ^j^^ 
variations for the flaring nights in more detail. An exponential function, e^^ (dashed line) is 
used to describe the rise and fall time for flux variations giving Xrise = 0.86 (x^/dof = 3.7/2) 
and Xfaii = —0.47 (x^/dof = 7.2/6), respectively. The characteristic flux doubling time, t^ar 
is estimated to be on the order of ~ 1 day. 
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Table 1. Summary of observations and analysis of lES 1218+304' 





Live Time 


Zenith 


Significance 


<!>(> 200 GeV) 


Units of Crab Nebula 




[hours] 


n 


['^] 


[10"" cm-2 s-i] 


flux (E > 200 GeV) 


2006-2007'^ 


17.4 


2-35 


10.4 


12.2 ±2.6 Stat 


0.05 ± 0.011 


2008-2009 


27.2 


2-30 


21.8 


18.4 ±0.9 Stat 


0.07 ±0.004 



Note. — Also shown arc the integral flux above 200 GeV assuming a spectral index F = 3.07, and 
the corresponding flux in units of integrated Crab Nebula flux over 200 GeV. 

''Including total live time, range of zenith angles for observations, and significance of the excess 

''(see 



Acciari et al. 



2009a^ 



Table 2. High-state and Low-state Data' 





Live Time 


/o 


r 




[minute] 






High 


195 


19.4 ±2.3 


3.01 ±0.17 


Low 


1436 


10.0 ±0.7 


3.12 ±0.10 


Avg 


1631 


11.5 ± 0.7 


3.07 ±0.09 



''The high-state and low-state data (sec the 
text for complete description) are well described 
by a power law fit of the form = /o ■ 

( o.5 Tev )~'"[ cm^°a Tev l- ""^'^ changes tO 

the VHE photon index during the fiaring activity 
are statistically insignificant. 



,xlO 



12 
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